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We present detailed chemical abundances for >20 elements in ^30 globular clusters in M31. 
These results have been obtained using high resolution (A/AA ^24,000) spectra of their inte- 
grated light and analyzed using our original method. The globular clusters have galactocentric 
radii between 2.5 kpc and 117 kpc, and therefore provide abundance patterns for different phases 
of galaxy formation recorded in the inner and outer halo of M3 1 . We find that the clusters in 
our survey have a range in metallicity of —2.2 <[Fe/H]< —0.11. The inner halo clusters cover 
this full range, while the outer halo globular clusters at R>20 kpc have a small range in abun- 
dance of [Fe/H]= — 1.6±0.10. We also measure abundances of alpha, r- and s-process elements. 
These results constitute the first abundance pattern constraints for old populations in M3 1 that are 
comparable to those known for the Milky Way halo. 



XII International Symposium on Nuclei in the Cosmos, 
August 5-12, 2012 
Cairns, Australia 



* Speaker. 

' The data presented herein were obtained at the W. M. Keck Observatory, which is operated as a scientific partner- 
ship among the California Institute of Technology, the University of California and the National Aeronautics and Space 
Administration. The Observatory was made possible by the generous financial support of the W. M. Keck Foundation. 



© Copyright owned by the author(s) under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike Licence. 



http: //pos .sissa.it/ 



The Detailed Chemical Abundance Patterns ofM31 Globular Clusters 



Janet E. Colucci 



1. Introduction 

Study of M31, the Milky Way's nearest massive neighbor, is interesting for many reasons. 
One of the most fundamental questions is whether the characteristics of the M31 spiral galaxy 
support the assertion that the Milky Way is a "normal" spiral galaxy. This is important because 
we can study the properties of the Milky Way in great detail, and studies of our own galaxy by 
necessity are the foundation for our understanding of how galaxies in general form and evolve. As 
the next closest massive galaxy, M3 1 is the first place to test galaxy formation theories developed 
from studies of the Milky Way, and in some respects is a more ideal test-case because M3 1 can 
be observed as a whole from the outside, whereas study of our own galaxy is complicated by our 
position within it. However, the distance to M3 1 means that we are unable to study its individual 
stars at the same level of detail that we can obtain in the Milky Way, using high resolution high 
signal-to-noise ratio (SNR) spectra and detailed chemical abundance analyses. 

With the development of our original technique for abundance analysis of high resolution 
integrated light (IL) spectra of globular clusters (GCs), we can now make significant advances in 
chemical evolution studies of distant massive galaxies. Unresolved GCs, which are luminous and 
therefore observationally accessible to large distances, can be used to learn about the chemical 
enrichment and formation history of other galaxies, just as they were originally used to learn about 
the formation of the Milky Way [|l|]. Our technique has been developed and demonstrated on 
resolved GCs in the Milky Way and Large Magellanic Cloud (LMC) in a series of papers ||, |], 

01. These works demonstrate that the IL analysis provides accurate Fe abundances and [X/Fe] 
ratios to ~0.1 dex, as well as distinguishes ages for GCs with a range in properties, including 
[Fe/H] of -2 to +0 and ages from 0.05 to 12 Gyr. 

With this method, we have now begun an unprecedented study of the chemical composition of 
the GC system of M31. Detailed abundances of ~ 20 elements were presented for a pilot sample 
of 5 M31 GCs in [Q]. As part of this ongoing project, here we extend the sample of [Q] and present 
ages, Fe, a-element and neutron capture abundances of an additional 22 GCs in M31. 

2. Observations 

We obtained high resolution IL spectra of the M31 GCs using the HIRES spectrograph on the 
Keck I telescope. The data were taken over several observing runs from 2006-201 1. In all observing 
runs we used identical setups that utilized the D3 decker, which has a slit size of 1.7" x 7.0" and 
spectral resolution of R = 24,000, which is sufficient to resolve individual spectral lines of GCs 
with velocity dispersions, a v >~ 710ns" 1 . The wavelength coverage is approximately 3800—8300 
A. Exposure times were between 1—5 hours for each GC. The SNR at 6000 A is approximately 
60-90. Data were reduced with standard routines in the HIRedux pipeline. 1 

3. Analysis 

As in all of our our previous spectroscopic analyses we first measure individual absorption line 
equivalent widths (EWs) using the semi-automated program GETJOB [g]. Line lists and oscillator 
strengths were taken from [pL ^J, ^] and references therein. In this work, we also use our IL spectral 

1 http://www.ucolick.org/ xavier/IDL/index.html 
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Figure 1: Left: The age-metallicity relationship. Black diamonds are from [|4p, and black squares are from 
this work. Right: The behavior of [Fe/H] with projected galactocentric radius from M31, where radii are 



taken from [ 1 1 1. Gray points are from [ 12 1, and measured using Lick system line indexes. 



synthesis code [J7|] for the best possible accuracy in our final measurements for Ca, Si, Eu, Ba and 
Y. Abundances for Ba and Eu include corrections for hyperfine splitting, as described in [^]. 

The IL abundance analysis is described in detail in []3|, |], |5]|. In this analysis we create syn- 
thetic GC color magnitude diagrams (CMDs) using the Teramo isochrones |Q] without convective 
overshooting, with extended asymptotic giant branch (AGB), and mass— loss parameter of 77=0.2. 
More discussion about the testing and choice of appropriate isochrones can be found in [Q, ||, [6||. 
Flux- weighted synthesized EWs of lines are calculated using our routine ILABUNDS [§J, which 



utilizes spectral synthesis routines from the 2010 version of MOOG [jlOJ.We use the ODFNEW 
and AODFNEW model stellar atmospheres from Kurucz. 2 All abundances are calculated under the 
assumption of local thermodynamic equilibrium (LTE). 

The age and [Fe/H] solutions for each cluster are identified as the range in synthetic CMD 
ages and [Fe/H] that produce the most self-consistent results using the 10—80 individual Fe I lines 
measured in each cluster. The best solutions have the smallest statistical error ((7n), and minimal 
dependence of Fe I abundance with line excitation potential (EP), wavelength, and EW. For each 
cluster there is a range in CMD ages that produce similarly self-consistent solutions. For older 
clusters this range is typically 10—15 Gyr, and leads to an internal systematic uncertainty in [Fe/H] 
of < 0.05 dex. For the total uncertainty in [Fe/H] for each cluster, we add this internal systematic 
age uncertainty, aA ge , in quadrature with the statistical error in the mean abundance, cTn/VN— 1, 
where N is the number of lines measured for each element. The best CMD solution is used to 
calculate abundances for the other elements. 



4. Results and Conclusions 

In Figure [I] we show the age-metallicity relationship for the 22 M3 1 GCs analyzed in this work 
and the 5 GCs analyzed previously [Q]. Figure [j] also shows the M31 GC [Fe/H] results against 
M31 projected galactocentric radii from [11]. The low resolution results of [12] are also shown 
for comparison and to highlight the greater precision of our technique when compared to low 
resolution results. We find a spread in [Fe/H] for GCs within Rm3i~20 kpc, and a fairly constant 
[Fe/H] for GCs at Rm31 >20 kpc. The mean value for the 5 GCs in our sample that lie at Rm31 >20 
kpc is [Fe/H]= —1.63 ±0.10. We therefore confirm the previous results of [13, 14], who also 
found a nearly constant [Fe/H] for outer halo GCs. Like our analysis, [ |14| ] found a flat metallicity 



http://kurucz.harvard.edu/grids.html 
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Figure 2: a-element abundances in our sample to date. Green squares show M31 GCs, Grey squares show 
MW IL abundances from our training set GCs [S], Small gray circles and orange triangles, respectively, 
show MW field stars [Il6|l and MW GC stars ifTvIl . 



1.5 
1.0 

0.5 

0.0 

-0.5 

-1.0 
-1.5 





i-B-i 




A 


Aa 

A 







-2.5 -2.0 -1.5 -1.0 -0.5 0.0 
[Fe/H] 




-2.5 -2.0 -1.5 -1.0 -0.5 0.0 
[Fe/H] 



Figure 3: [Ba/Fe] and [Y/Fe], symbols are the same as in Figure 2. 
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Figure 4: [Eu/Fe] and [Eu/a], symbols are the same as in Figure 2. 



distribution for GCs outside 30 kpc, with a mean metallicity of [Fe/H]=— 1.6, which was estimated 
from line indexes. We also note that at least one GC in our sample has an intermediate age (1-3 
Gyr) at high [Fe/H]; this cluster has disk-like kinematics [15] and lies within 7 kpc. 

In Figures ||-[|, we show the [X/Fe] as a function of [Fe/H] for Ca, Si, Ba, Y, and Eu. For 
the a-elements Ca and Si, we find a plateau value that is similar to the value in MW GCs and 
field stars from flip], [l7|], which is consistent with rapid early star formation in M31. Note that 



our results do not confirm the line index results from [fi~8l Il9|]. The underestimations of [a/Fe] in 



[18, 19] are likely due to Mg depletion within some of the cluster stars; we first published evidence 
of this effect being visible in the integrated light of GCs in [|], ^]. The star-to-star variations in the 
light elements Na, O, Mg, and Al, are a well known phenomenon, and discussed in the review by 
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Figure 5: [Ba/Y] and [Y/Eu], symbols are the same as in Figure 2. 



[20]. Most recently [ pi| ] confirmed, using a slightly different IL analysis, that star-to-star variations 
in Mg are detectable in IL, and therefore Mg is not a good proxy for [a/Fe] in IL GC spectra. In 
Figure || there is a clear knee visible in the [Ca/Fe] values, which closely matches the abundance 
patterns of MW field stars. Obvious outliers in the abundances of our GC sample include at least 
one GC with high [Fe/H] and high [a/Fe], and one GC with low [Fe/H] and low [a/Fe], the latter 



of which may be associated with an accretion event [22]. 

Note that the abundances shown in Figures 3 and 4 are the first results for r- and s-process 
abundances in old populations in M31. In general, we find that the Ba, Y, and Eu are consistent 
with MW heavy element abundances. The [Eu/Fe] values for GCs with [Fe/H]< —1.5 show some 
suggestion of being higher than the mean of MW GCs, although they are still within the range 
exhibited by the MW. A larger sample of low metallicity Eu abundances is necessary to confirm 
this result. 

In Figure |-|we show key neutron capture abundance ratios as a function of overall metallic- 
ity. These are particularly important diagnostics as [Ba/Y] identifies the yield of second peak to first 
peak s-process elements, while [Y/Eu] identifies the yield of light s-process to r-process elements. 
[Eu/a] is a diagnostic for different elements that are primarily produced in Type II supernovae. We 
find that this sample of M3 1 GCs is dominated by r-process enrichment at low metallicity, and the 
s-process influence appears to increase around [Fe/H]~ —1, similar to the MW. These results, like 
the [a/Fe] results, are consistent with rapid, early star formation in M31. 

In summary, we present ages and abundances of Fe, Ca, Si, Ba, Y, and Eu for a total sample of 
27 GCs in M31. For the abundances shown here, M31 generally looks like the MW; even this basic 
similarity has been controversial to assert from previous line index analyses. Our analysis of the 
IL spectra of extragalactic GCs demonstrates the importance of quantitative, accurate abundance 
ratios in clarifying the chemical enrichment and formation history of galaxies beyond the Milky 
Way. M31 is the first of many galaxies within the Local Group and nearby Groups that we are 
studying with this technique. 
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